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Abstract 

We present a method for the analytic solution of small x structure functions. The 
essential small x logarithms are summed to all orders in the anomalous dimensions and 
coefficient functions. Although we work at leading logarithmic accuracy, the method is 
general enough to allow the systematic inclusion of sub-leading logarithms. Results and 
predictions are presented for the gluon density, and the structure functions F2(x,Q 2 ) 
and Fl(x,Q 2 ). We find that corrections to the simple double logarithmic calculation are 
important in the HERA range and obtain good fits to all available data. 



As a result of the recent work of Catani and Hautmann |y, it is now possible to include 
the dominant small x dynamics encompassed by the formalism of Balitsky, Fadin, Kuraev and 
Lipatov (BFKL) [0 within the framework of the renormalisation group and collinear factorisa- 
tion, and some (mostly numerical) studies have already been performed ||. In this paper, we 
wish to present an analytic solution to the relevant evolution equations and their convolution 
with the appropriate coefficient functions. Throughout we work in the high energy limit, i.e. 
we sum all terms in the perturbative expansion of the cross section which are 




where s is the relevant centre-of-mass energy and Q 2 characterises the typical short distances 
involved. We shall focus on deep inelastic scattering at the DESY ep collider, HERA. In which 
case, \fs is the 7p centre-of-mass energy and — Q 2 is the photon virtuality, i.e. the Bjorken-x 
variable, x ~ Q 2 / s. Our approach is quite general and it will be clear how to extend it beyond 
the leading logarithmic accuracy. 



Altarelli-Parisi Evolution at small x and the gluon density 

Recall the Dokshitzer, Gribov, Lipatov, Altarelli, Parisi (DGLAP) equations for the parton 
distribution functions ||: 

d -^^ = T,iAUQ 2 ). (i) 

fhiQ 2 ) i s the Nth moment of the momentum distribution for partons of type i and Yn * s the 
anomalous dimension matrix, i.e. 



Pn(Q 2 ) = f'dxx^f&Q 2 ), 
Jo 

7jv = / dxx N Pij(x). 

•/ 



Our notation is such that the important gluon anomalous dimension, 



a, fa, N 1 



7^f = ^ + 2C(3)^j +.... (3) 
and a s = 3a s /7r. 

These equations are solved given the boundary conditions, f^(Qo), i.e. they allow the 
<5 2 -dependence of the parton distribution functions to be determined but not their absolute 
normalisation. 
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In the high energy (i.e. small x) limit, we keep only those terms in the anomalous dimension 
matrix which are ~ (a s /N) n , i.e. the leading logarithmic terms in the splitting functions. In 
this case, evolution is driven by jff , which satisfies 0, § 

1 = jfxitf), (4) 

where 

x(7) = 2V(l)--0(7)--0(l-7) (5) 
and if>('j) * s the Euler-gamma function. The first two non-zero terms in the series expansion 
are written in eq.(3). The DGLAP equations then have the simple solution: 

f s N (Q 2 ) = fKQl), 

- - --- 4 -^e X p(r^)-^ (6) 



fUQ 2 



fUQl) + ~f s N(Ql) 



The singlet quark density is f$j{Q 2 ) = J^ifhiQ 2 ) where the sum runs over all quarks and 
anti-quarks. Since we work in the small x region, we expect the gluon density to be dominant 
and subsequently drop all reference to the singlet density (except implicitly in the input to 
F 2 (x, Q 2 )). We have explicitly checked that this makes very little quantitative difference to our 
results. 

In order to construct a sensible gluon structure function, we do not merely invert the N- 
space solution above. It is more natural to define the gluon structure function to be that object 
which would be observed if we had a coloured current available as our probe. In which case 
there are important contributions which arise, not only from the QCD evolution but also from 
the coefficient function. One can think of such corrections as arising from graphs which should 
not be exponentiated via the renormalisation group and so contain no explicit strong ordering 
of the rung momenta. These graphs are essential for a sensible definition of the gluon structure 
function (e.g. as the object which is closely related to the structure function Fl(x, Q 2 )) and for 
consistency with the structure function which is constructed by integrating the 'unintegrated 
gluon density', JF(x, k 2 ), obtained by solving the BFKL equation 0. To see this, we start from 
the BFKL definition of the gluon structure function, i.e. 

G(x, Q 2 ) = f° ^ f 1 -5(1 - x/z)Q(Q 2 - k 2 )T(z, k 2 ). (7) 

JO rC Jx Z 

The hard scatter cross section which is to be convoluted with F(x, k 2 ) is thus 

v N (k 2 /Q 2 ) = f 1 dxx N -H{l - x)Q(Q 2 - k 2 ). (8) 
Jo 



Now, in the formalism of Catani and Hautmann [|TJ, 



r°° dk 2 ( P \ 7jv 

G N (Q 2 ) = i 9 N 9 1 f q5 J Mk 2 /Q 2 ) Rn f N (Q 2 ) 



i.e. 



G N m = r n fUQ 2 ' 



(9) 



(10) 



Here Rn plays the role of a coefficient function for the 'natural' definition of the gluon structure 
function |7], This important factor, which is a process independent but scheme dependent 
quantity is given, in the MS scheme (which we shall employ throughout this paper) by 

f r(i - -y#) x(t#) ) 1/2 



R 



N 



lr(i + 7#) [-7^fx'(7f)] 
exp {nffiil>(l) 



X 



7 - ^ y( 7 )-^(i-7) 
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Solution in x-space 

Let us now show how to invert the solution for Gjq{Q 2 ) back into x-space. We must perform 
the integral 

G(x,Q 2 ) = -L /" dNf^(Ql) ^ exp(iVlnl/x + ^(g 2 

27TZ J v 

where the integral is over a contour to the right of all singularities, and 

rQ 2 dk 2 



'12) 



Z N (Q 2 ) 



99 



2 TN 



Ql k 

11+ y 



r 



' o> s (Q 



m—l 



N 



' a s (Q 2 
N 



m—l 



E 



MO 



i N r 



(13) 



( = \ia(a s (Ql) / a s (Q 2 )) , ■y 2 = 12//3 and, for consistency with the standard approach, we run 
the coupling at the scale k 2 (in the anomalous dimension integral). However, we note that at 
the leading logarithmic accuracy this is an essentially free choice. The coefficients a n define the 
series expansion of the gluon anomalous dimension: 



In 



oo 



n=l 



NJ 



(14) 
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We write the factor also as a series expansion: 



oo / 1 \ 

n=0 ViV 7 



(15) 



and we note that R N = 1 + §C(3)(a s /iV) 3 + .... 



We choose the boundary condition, 



f 9 N (x,Q 2 ) = AfQ(x -x), 



which becomes f^iQl) = N ' /N in moment space if we choose x/xo as the conjugate variable 
to N. This starting distribution is a good approximation to the expected ~ (1 — x) 5 behaviour 
if Xq ~ 0.1, and leads to reliable results for x < 0.01. The choice of a flat starting distribution 
(at small x) is motivated by the known behaviour of total cross sections at high energies, i.e. 
the 'soft' pomeron is known to have intercept close to 1 0. It is the small x behaviour one 
would expect in the absence of any perturbative QCD corrections. 

We can now perform the iV-plane integral by taking a particular choice of contour to be the 
line from r — zoo to r + zoo plus a circle with centre at the origin and radius r > 4a s (<5o) In 2. 
This choice of contour ensures the analyticity of the integrand along the contour, i.e. that all 
cuts lie within the circle (e.g. for the precise location of the cuts see ref.[[|). The value of the 
integral is now equal to that over the circle, and putting N = re ld we obtain in a straightforward 
manner: 



where £ = ln(xo/x). We note that exactly the same method could be used if we were to choose a 
powerlike input, or even the (1 — x) 5 behaviour. One simply finds the moment space expression 
for the input and expands in powers of N. We also note that for small x the result obtained using 
the saddle-point method to evaluate eq.(12) does not give a good approximation and provides 
misleading results. This failure occurs essentially because, along the contour of steepest decent, 
the integrand does not fall quickly enough for values of N far from the saddle-point. 
Let us now discuss our solution. Firstly, we see explicitly the double log result, 




i/2 



(16) 



G(x,Q 2 )~/ (2 7v /eC), 
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which arises when only the leading order (in a s ) terms are kept. Going beyond this first term, 
the inverse factorials associated with the Bessel functions ensure that the summations in our 
expression converge for all x, despite the fact that the expansions of 7^ and Rn diverge for 
iV < 4a s In 2. This effect of convergence in x-space was pointed out using a similar, but slightly 
less direct, argument in |J. Using the first sum in eq.(16) we can recover the power behaviour 
of the structure function at small enough x. It arises after taking the small argument expansion 
of the Bessel function (which is appropriate for large order Bessel functions) and using the fact 
that c n+ i/c n = 4a s ln2 for asymptotically large n, i.e. the general term in the sum over % is 



A similar power behaviour is generated by the other terms in eq.(16), i.e. due to the Q 2 - 
evolution, but they are not important for all practical values of x. The same qualitative 



approximate solution for f 9 (x, Q 2 ). In other words, the dominant corrections to the double log 
result are due to the presence of the Rn factor, i.e. the corrections to the evolution are small 
(due to the relatively small size of the coefficients in the expansion of the gluon anomalous 
dimension), only becoming dominant at very large Q 2 and/or very small x. 

The fraction of G(x, Q 2 ) which arises solely from the double log graphs (i.e. the Iq Bessel 
function) is presented in the contour plot shown in fig.(l). It can be seen that the high energy 
(BFKL) corrections are significant over the HERA range despite the fact that the coefficients c\ 
and C2 vanish. We note that the contribution from the BFKL corrections to the evolution (i.e. 
those terms involving the hi coefficients) are almost entirely negligible, in fact they contribute 
less than 4% over the x-Q 2 range probed at HERA. In fig. (2), we show the x dependence of 
G(x, Q 2 ) at different Q 2 values and compare to the double log contribution. In all of our plots, 
we choose xq = 0.1 and take M = 1.1 and Qq = 2.0. M and Qq are the only parameters for 
the gluon, and are fixed by fitting F 2 (x, Q 2 ) to the HERA data (see the following section for a 
discussion of this procedure). Note that our approach does not permit a flat gluon structure 
function, even though our input density was flat. This is in keeping with the standard BFKL 
result developed by direct solution of the BFKL equation. The scale Qq is to be understood 
as the scale below which we cannot use the perturbative approach. As such, we are unable to 
make any definite statements regarding the eventual saturation and flattening off of the small 
x structure functions (as Q 2 falls below Qq) since this procedure is governed by physics beyond 




conclusion has been reached in ref.[PL where a sum of Bessel functions was presented as an 
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that which is considered here. (However, we do see a hint of the breakdown of our approach, 
as we will discuss in the next section.) Indeed, we shall find that the quality of our fits to 
F 2 (x,Q 2 ) is largely insensitive to the choice of Ql once it is above ~ 1 GeV 2 . Our approach 
should be contrasted to that which attempts to evolve from some, typically quite low, value of 
Qq with a flat (or valencelike) starting distribution to higher Q 2 ]10|, [II]. 



Deep inelastic structure functions 

In the previous section we concentrated on the gluon structure function, G(x,Q 2 ). It involves 
no new techniques to extend the formalism to the case of the deep inelastic structure functions, 
F 2 (x, Q 2 ) and F L (x, Q 2 ). Eq.(lO), which defines the gluon structure function, is merely a specific 
form of the more general expression for the dimensionless cross section, F(x, Q 2 )'- 

F N {Q 2 ) = C N (Q 2 /f/ F )fM) (17) 

where ji F is the factorisation scale (chosen to equal Q 2 ) and Cn(Q 2 / is the coefficient 
function (equal to Rn in the case of the gluon structure function). 

Catani and Hautmann have shown that the coefficient function can be factorised into a 
product of the process independent (but factorisation scheme dependent) factor, R N , and a 
process dependent factor, /ijv(7|f)j where [Q 

h N (l)^l r^-(K-Y Mk 2 /Q 2 )- (18) 



lo k 2 \fi 2 F/ 

The hard subprocess cross section, a^{k 2 /Q 2 ), is the lowest order (in a s ) cross section for 
scattering off-shell gluons (off the virtual photon in the case of deep inelastic scattering). 
Thus, for the structure function, Fl(x,Q 2 ) 

F L , N (Q 2 ) = <& )h L Ai 9 N 9 ) Rn f N {Q 2 ) (19) 

with 

h M- as NT 4(1-7) r 3 (l- 7 )r 3 (l + 7 ) 

^(7) - ^N f T R 3 _ 27 r(2 _ 27)r(2 + 27) (20) 

and where (e 2 ) is the mean quark charge squared. So, in order to evaluate F L (x,Q 2 ), we 
merely replace the c n coefficients in eq.(16) by the corresponding coefficients in the expansion 

Of h LtN (lN) R N- 
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Similarly, 



h2,N{V = 2(1 _ 7 ) L ^ ^ ^ 



determines the Q 2 -dependence of F 2 (a;, Q 2 ), i.e. 



«91nQ 2 



(4)[^(a s (g 2 )) 7 f + 2n /7 M(g 2 ) 

(el) h, N (7 9 N 9 ) Rn f 9 N (Q% (22) 



to formally leading order. 

In fig. (3), using the same choice of parameters (no more are needed) as in the discus- 
sion of G(x, Q 2 ) previously, we present our predictions for the longitudinal structure function, 
Fl(x, Q 2 ). As well as the full solution, we show the double log contribution to it. The largeness 
of the corrections to the double log calculation (in comparison to case of G(x,Q 2 )), can be 
traced back to the fact that the second and third coefficients in the expansion of Hl^nRn are no 
longer zero. Also shown in fig. (3) is the result of re-fitting the HERA data on F 2 (x, Q 2 ), while 
keeping only the double log Bessel function. In order unambiguously to establish the existence 
of the high energy corrections, it is ultimately necessary to expose deviations from the double 
log approach (or more precisely approaches which do not sum the infinity of high order cor- 
rections 0(a s /N)) and so this is the reason for our comparison. As seen, the prediction from 
the double log approach is mostly larger than that for the full expression, but flatter with x. 
This largeness comes about mainly because the starting scale is much lower, and hence there 
has been more time for evolution to take place. 

Let us now turn to the structure function, F 2 (x, Q 2 ) and its comparison with the HERA 
data [12|. We start by considering the expression of eq.(22). In order to construct F 2 (x,Q 2 ), 
we must integrate over Q 2 and invoke an input distribution, F 2 (x,Qq). We choose this to be 
of the form A + Bx~ x . We see no reason to believe that the input form of F 2 (x, Q 2 ) should 
be purely flat since, as demonstrated, the gluon structure function always has some powerlike 
behaviour due to the coefficient function. Indeed, we are not able to obtain a very good fit with 
a completely flat input. We could have chosen an input of Bx~ x , and still obtain a comparable 
fit. However, our aim is not simply to obtain the best fit with the least number of parameters, 
but to determine the behaviour of the structure function as accurately as possible, and we 
believe the chosen input is the best way to do this. This introduces three extra parameters. 
The values of our 5 parameters (Qq, the normalisation of the input gluon density and the three 
parameters in F 2 (x,Qq)) are then obtained by fitting to the HERA data. Throughout, we 



work with 4 quark flavours and Aqcd = 115 MeV. Agreement with the data is very good, i.e. 
X 2 = 48 for the 92 data points which have x < 0.01. 

Although we obtain very good agreement with the available HERA data, we expect our 
results to be subject to important corrections. Let us now explain why. In the leading log 
approximation, the structure function F 2 (x, Q 2 ) can be written in the form 

F2MQ 2 ) = (e^iMQ 2 ))/^ ^ 

+ ln f {e 2 q )C% f %{l q ^W) + i$fM)- (23) 
J Qo 1 

The superscript on the coefficient functions specifies the order (in a s ) of the contribution, i.e. 
Ci.at = SnCi,7v- The parton model coefficient functions are C^ N = and C^ N = 1. 

Taking the derivative of this expression leads to eq.(22), but only after neglecting the higher 
order terms which are induced by differentiating the coefficient functions. Such terms are 
formally sub-leading since 

d 2 f3 d 



-a 



dhiQ 2 s 4tt da s ' 

However, they are not sub-leading once eq.(22) is integrated to form the structure function, 
F 2 (x,Q 2 ). 

To see how important these corrections are expected to be, we expand the coefficient function 

<$r=E^ «.(!)"■ (24) 

n=0 ViV 7 



The ratio of the term ~ a s (a s /N) m in dC g N / \n.Q to the corresponding term in the series 
expansion of C^ N 7^f is thus 

/ (3 a s \ (m + l)b m ro ., 
— TT — " • ( 25 ) 

Since a n+ i/a n = 4 In 2 for large n and, assuming a similar relation for the p n coefficients, it 
follows that this ratio becomes 

(3 a s 1 

rsj 

47T A 

where a n m A (41n2) n_1 . Since A < 1 we cannot ignore such contributions. We should 
emphasise that p n +i/p n cannot exceed 4 In 2 (since we know the dominant singularity arises at 
N = 4a s ln2) and that assuming any p n +i/p n < 4 In 2 leads to an even stronger enhancement 
of the derivative terms (e.g. by a factor of m for p n+ i/p n 4 In 2). All the evidence from the 
calculation of the series expansion of the coefficient function is that p n +i/p n is indeed ~ 4 In 2 
for large n. 
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Although these corrections are formally sub-leading, we believe it is unlikely that they will 
cancel with higher order graphs and as such it is safe (and probably more appropriate) to 
include them in this order of the calculation. Since they lead to a negative contribution to 
dF 2 (x, Q 2 )/dliaQ 2 , we therefore expect a reduction in dF 2 (x, Q 2 )/dhiQ 2 as x falls at fixed Q 2 
(or Q 2 falls at fixed and small x). 

Using the expressions of ref . we have computed the series expansions of the gluon coeffi- 
cient function, C^ N , and the anomalous dimension, 7^, in the MS scheme to 18th order (we 
will shortly have more to say on the choice of scheme), the first six terms of course agreeing 
with the numerical values explicitly given in ref.0. It is then a simple matter to calculate F 2 by 
transforming eq.(23) (using eq.(6) to give f^(Q 2 )) to x-space, in the same way that we obtained 
eq.(16) from eq.(12). We make the simplification that (as in eq.(16), and as when obtaining 
Fl(x,Q 2 ) by transforming eq.(19) to x-space) we drop the quark singlet density in the defini- 
tion of the gluon distribution (as explained in the remarks following eq.(6)). We also neglect 
7jv and the 0(a s ) contributions to the quark coefficient function when calculating F 2 (x, Q 2 ) 
(since the input quark density is small in comparison to the evolved gluon density, the quark 
coefficient function is smaller than that for the gluon and 7^ is smaller than 7^). Using this 
method, we compute F 2 (x,Q 2 ) including those corrections which were neglected when eq.(22) 
was integrated over Q 2 . The solid line in fig. (4) shows the result of a new fit to the HERA data 
and marginal improvement in the \ 2 °f the fit is found, i.e. \ 2 = 45 for the 92 data points. 
A considerable improvement in the insensitivity to the value of Qq is also found. For our best 
fit Qq = 2.0 and M = 1.1, and these are the parameters used to determine our predictions for 
G(x, Q 2 ) and F L (x, Q 2 ). We also find the input to F 2 {x, Q 2 ) to be 0.15 + 0.035x" a4 . 

From our results we conclude that our choice of a G-function form of the gluon input is 
appropriate. Also, we note that the x 2 onr y starts to worsen significantly once Ql < 1 GeV 2 . 
This is consistent with idea that the scale at which we define our input should be essentially 
arbitrary, providing it is large enough for the perturbative expansion to apply, and not too large 
to fill the available phase space. The dotted line in the figure shows the previously discussed 
best fit for F 2 (x, Q 2 ), i.e. ignoring the derivatives of the coefficient function. The dashed 
line shows the result of fitting the data taking only the leading term in the Bessel function 
expansion, i.e. the double log result, and flat inputs for the gluon and for F 2 (x, Q 2 ). This also 
has a very good \ 2 °f 44 for the 92 points, but high sensitivity to Q\. It is clear from the plots 
that the differences between the full leading log calculation of F 2 (x, Q 2 ) and the dotted (and 
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dashed) line are consistent with our expectations. 

As a slight word of caution we must note that for low enough x, dF 2 (x, Q 2 )/dlnQ 2 even- 
tually becomes negative, and F2(x,Q 2 ) rises for falling Q 2 . It is this region of negative \nQ 2 
derivative where our calculation starts to become untrustworthy. However, the effect only sets 
in for x ~ 10~ 4 at Q 2 = Qq, and the value of x at which dF 2 (x, Q 2 )/dln Q 2 = falls very rapidly 
as Q 2 increases. Similarly, we feel that the predictions for the gluon and Fl(x,Q 2 ) should be 
viewed with a little caution at extremely low x and small Q 2 . Nevertheless, it is reassuring 
that the region of breakdown is where we might expect physics beyond that considered in our 
approach to become important. 

It is also important to comment on our choice of scheme. We could just as well have 
computed F 2 (x, Q 2 ) in the DIS scheme, and obtained precisely the same results as in MS (see 
below). This is true providing we take care to include the sub-leading corrections to 7^ which 
contribute in the leading order to F 2 (x,Q 2 ). These terms are those neglected in eq.(5.27) of 
Catani and Hautmann when transforming to DIS scheme. 

Before concluding, we wish to make a few remarks regarding the recent conclusions of Ellis 
et al [f|. Recall that F 2 ^ N {Q 2 ) — F 2N {Q^) ~ a s O(a s /N), and as such any calculation which 
sums the leading logarithms must necessarily include the quark anomalous dimensions, 7^ 
(and 7jv), computed to the same accuracy (i.e. ~ a s O(a s /N)). Thus, it is not sensible to 
compute F 2j n(Q 2 ) by including only the BFKL anomalous dimension, 7^ ~ 0(a s /N), and 
(for example) the two-loop form for the rest of the anomalous dimension matrix. It is also 
important to appreciate that, in the DIS scheme (where all the gluon coefficient functions are 
zero) the quark anomalous dimension 7^ is where the physics which gave the dominant (leading 
log) contribution to the gluon structure function G(x, Q 2 ) resides, i.e. 

2n f 7$ = R N h 2 ,M 9 ) (26) 

determines the DIS scheme anomalous dimension in terms of Rn (calculated in MS scheme). 
This is in contrast to the MS scheme, where it is the coefficient function C^ N which has the 
large coefficients in the series expansion in a s /N (due to the corresponding large coefficients in 
Rn) and, as a result, the coefficients of 7^ are much smaller in the MS scheme than in the DIS 
scheme. Therefore, the relatively small effect observed in ref.|l when only the BFKL anoma- 
lous dimension is added to the full two-loop formalism is consistent with our finding that the 
bi coefficients in eq.(16) have little impact. However, on adding the all orders quark anomalous 



10 



dimension, 7^, Ellis et al found huge effects which are very sensitive to how momentum conser- 
vation is applied. The largeness of their effect is driven by the largeness of the coefficients in the 
expansion of 7^ as calculated in the DIS scheme. We wish to urge caution in interpreting their 
results as evidence of important higher order corrections to the leading log BFKL framework. 
In particular, as we have just emphasised the bulk of the BFKL physics resides in the leading 
log factor R N (i.e. in the DIS scheme 7^), so part of the large effect which is seen after turning 
on 7jy is actually a leading effect. In addition, on adding the quark anomalous dimension the 
leading order expression or F 2 (x, Q 2 ) is given by eq.(23), whilst the next-to-leading order result 
is more complicated and involves, in particular, terms proportional to (7^f) 2 . Such terms are 
small in the MS scheme but may be large in DIS. The corresponding large terms never appear 
in a MS scheme calculation since the coefficient function, C^ N , is never iterated (i.e. raised to 
a power). Consistency between the schemes therefore suggests a cancellation of large terms in 
the DIS scheme. In order to perform a complete next to leading order calculation of F 2 (x, Q 2 ) 
there are many other terms which need to be computed, i.e. the ~ a s O(a s /N) corrections to 
the gluon anomalous dimensions 7^ (and 7^) and the higher order corrections to the quark 
anomalous dimensions, i.e. ~ a 2 0(a s /N). Cancellations should then occur between the terms 
proportional to the square of the leading order quark anomalous dimension and these other 
terms. By not performing the full next-to-leading calculation, Ellis et al are unable to observe 
these expected cancellations and as such their conclusions, based upon different assumptions 
regarding the imposition of momentum conservation, may be drawn prematurely. 

Conclusions 

We have presented an analytic approach to the evaluation of small x cross sections and studied 
the behaviour of the gluon structure function, defined in a way which is consistent with the 
previous studies based upon the direct solution to the BFKL equation. In particular the solution 
in x-space is obtained exactly, and explicitly reveals the extreme limitations of the double 
leading log result and the power behaviour expected from the BFKL approach. In addition, we 
examined the deep inelastic structure functions F 2 (x, Q 2 ) and F L (x, Q 2 ) and demonstrated that 
the high energy corrections (to the double log calculation) are significant in the HERA region. 
Consistency with the data on F 2 (x, Q 2 ) is found. The structure function F L (x, Q 2 ) and the Q 2 
dependence of F 2 (x, Q 2 ) should be able to provide sensitive tests of the small x dynamics; in 



II 



particular deviations from the traditional approach (expansion in a s ) may well be observable. 

We have not discussed the process of heavy quark production (in deep inelastic scattering or 
in photoproduction), although this process also ought to shed important light on the essential 
dynamics [§. Also, the recent measurement by the ZEUS collaboration of the dijet cross section 
in photoproduction [13[] could be confronted with theory using the techniques presented here 
]T4j. Finally, we wish to make available the expansion (in a s /N) of the coefficient function, R^, 
and the MS scheme expansions of the quark anomalous dimension 7^ and the gluon coefficient 
function C^ N . These are displayed in the following table. 
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Table 



Values of the first 18 coefficients in series expansions: 



<J<?(MS) _ a s 1 ^ q 

lN ~27c3to an V'.vV • 
^(i)(Ms) ai s 1 ^ _ , , , ^ 



n=0 



oo 

>(MS) 



71=0 ViV 7 

i-e- Pn = yfPn (see eq.(24)) and c„a™ = c n (see eq.(15)). 
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a q9 


T) 

fn 


c 


o 


1.00 


1.00 


1.00 


1 


1.67 


1.49 


0.00 


2 


1.56 


9.71 


0.00 


3 


3.42 


1.64X10 1 


3.21 


4 


5.51 


3.91X10 1 


-0.811 


5 


7.88 


1.29xl0 2 


4.56 


6 


2.57X10 1 


2.41xl0 2 


3.27X10 1 


7 


4.42X10 1 


6.53xl0 2 


-2.95 


8 


8.77X10 1 


1.93xl0 3 


1.08xl0 2 


9 


2.83xl0 2 


4.01xl0 3 


4.00xl0 2 


10 


5.11xl0 2 


1.14xl0 4 


1.33xl0 2 


11 


1.24xl0 3 


3.17xl0 4 


2.10xl0 3 


12 


3.68xl0 3 
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Figure Captions 

[1] Contour plot exhibiting the contribution to the full gluon structure function made by the 
double leading log term. 

[2] Gluon structure function G(x, Q 2 ) as a function of x plotted for a range of Q 2 values. 
Contribution made by double leading log approximation shown by dot-dashed line. 

[3] Prediction of longitudinal structure function Fl(x,Q 2 ) as a function of x plotted for a 
variety of Q 2 values. Contribution made by double leading log term shown by dot-dashed 
line. Also shown by dashed line is the prediction made using the best fit for F 2 (x,Q 2 ) 
while keeping only the double log term. 

[4] Comparison of theoretical predictions with the small x (i.e. x < 0.005) data from the 
(a) ZEUS collaboration (renormalised up 2%) and the (b) HI collaboration (renormalised 
down 4%). The dotted line corresponds to the best fit for this expression minus the 
formally subleading terms coming from the derivative of the coefficient function and the 
dashed line to the best fit for the double leading log approximation. 
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